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We report the first measurement of strange (A) and anti-strange (A) baryon production from 
^s^j^ = 130 GeV Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC). Rapidity density 
and transverse mass distributions at mid-rapidity are presented as a function of centrality. The 
yield of A and A hyperons is found to be approximately proportional to the number of negative 
hadrons. The production of A hyperons relative to negative hadrons increases very rapidly with 
transverse momentum. The magnitude of the increase cannot be described by existing hadronic 
string fragmentation models. 

PACS numbers: 25.75.Dw 



Ultra-relativistic nucleus-nucleus collisions provide a 
unique means to create nuclear matter of high energy 
density (temperature) and/or baryon density over an ex- 
tended volume . The first results from the Relativistic 
Heavy Ion Collider (RHIC) have shown that the large 
charged particle multiplicity, measured in Au-|-Au colli- 
sions at y/s^ — 130 GeV, corresponds to an energy den- 
sity significantly higher than that previously achieved in 
heavy ion collisions |[ In addition, the anti-proton 
to proton ratio at mid-rapidity has been measured to be 
in the range 0.6-0.7 [|[ ^, 0, which is indicative of par- 
ticle production from a low net baryon density regime. 
Thus the global characteristics of nucleus-nucleus colli- 
sions at RHIC are the formation of a high energy density 
and low net baryon density region at mid-rapidity. 

The yield of baryons and anti-baryons produced in rel- 
ativistic nuclear collisions is very important because it is 
sensitive to two fundamental, not yet fully understood 
aspects of hadron production dynamics: baryon/anti- 
baryon pair production and the transport of baryon num- 
ber from beam rapidity to mid-rapidity. The nature 
of baryon/anti-baryon production itself is the subject 
of much interest, with theoretical conjecture address- 
ing a range of possible mechanisms from string frag- 
mentation B, to exotic mechanisms involving Quantum 



Chromo Dynamics (QCD) domain walls The phys- 
ical nature of the entity which carries baryon number, 
and the means by which baryon number is transported 
over a large rapidity gap into the mid-rapidity region are 
also subjects of considerable experimental and theoreti- 
cal interest 0, iii|, m. 

Strange baryon/anti-baryon production is particularly 
interesting due to the increased sensitivity to the avail- 
ability of strange/anti-strange quarks, which is expected 
to be suppressed relative to light quarks in hadronic mat- 
ter due to the strange quark mass. Strangeness produc- 
tion has long been predicted to be a signature of Quark 
Gluon Plasma (QGP) formation ||l^. The strangeness 
production in previous generations of heavy ion exper- 
iments has been observed to be significantly increased 
compared to those from p-fp, p-f A and light ion colli- 
sions [|lj, |l^, |l^, although questions remain about 
the exact strangeness production mechanism. In partic- 
ular, the relative importance of strange baryon produc- 
tion from hadronic rescattering differs between calcula- 
tions |l8[ |l^, depending on both the evolution of the 
system and the scattering cross sections assumed. Ex- 
otic dynamical mechanisms that have been proposed for 
strange baryon production include, for example. Color 
String Ropes ||2^, String Fusion and Multi-mcsonic 
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Reactions ||2^. All require a high local energy density 
and therefore suggest that strangeness production occurs 
early in the collision. 

In this letter we report on mid-rapidity (|?/|<0.5) 
Lambda (A) and anti-Lambda (A) production in Au-|-Au 
collisions at ^s„„ — 130 GeV. The data were taken 
with the STAR (Solcnoidal Tracker At RHIC) detec- 
tor. The main components of the detector system for 
this analysis have been described in detail elsewhere . 
They included a large volume Time-Projection Chamber 
(TPC) H, a pair of Zero-Degree Calorimeters (ZDC) lo- 
cated at +/- 18 meters from the center of the TPC, and a 
Central Trigger Barrel (CTB) constructed of scintillator 
paddles surrounding the TPC. 

The TPC was used to provide tracking information 
and particle identification by measuring the specific ion- 
ization loss (dE/dx) for charged particles. The ZDCs 
were used in coincidence to define a minimum bias trig- 
ger by selecting interactions from the intersection region 
of the colliding beams, and to measure the forward go- 
ing energy of near-beam rapidity fragments, consisting 
primarily of spectator neutrons. The CTB was used to 
define a central interaction trigger based on the energy 
deposition from charged particles entering its acceptance. 

Data from both the minimum bias trigger and central 
trigger have been used for this analysis. Similar to a pre- 
vious analysis 0| , the collision centrality was defined off- 
line using the total charged particle multiplicity within a 
pseudo-rapidity window of \r]\ < 0.5. The charged parti- 
cle multiplicity distribution was divided into five central- 
ity bins, corresponding to approximately the most cen- 
tral 5%, 5-10%, 10-20%, 20-35%, 35-75% of the total 
hadronic inelastic cross section of Au+Au collisions. 

The A and A particles were reconstructed from their 
weak decay topology, A ^ p7r~ and A pi^^ , using 
charged tracks measured in the TPC. These tracks were 
projected back from the TPC to find their Distance of 
Closest Approach (DCA) to the primary vertex. Since 
the tracks of particles from A and A decays should not 
originate from the primary vertex, only proton (anti- 
proton) candidate tracks missing the primary vertex by 
0.90 cm and tt" (7r+) candidate tracks missing the pri- 
mary vertex by 2.85 cm were selected. Particle assign- 
ments for p (p) and 7r~ (tt"'") candidates were based on 
charge sign and the mean energy loss, <dE/dx>, mea- 
sured for each track. Candidate tracks were then paired 
to form neutral decay vertices, which were required to be 
at least 5 cm from the primary vertex. The momentum 
vector of the reconstructed parent particle was required 
to originate from the primary interaction: the distance of 
closest approach of the parent trajectory to the primary 
vertex was required to be less than 0.5 cm. The effect 
of these cuts on the reconstruction efficiency was deter- 
mined using the embedding techniques described below. 

Figure |l| shows the invariant mass distributions for the 
reconstructed A and A candidates in < 0.5 for two typ- 
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FIG. 1: Invariant mass distribution of A (solid-circles) and 
A (open-circles) candidates from the data sample analyzed 
including both central and minimum-bias collisions: (a) 0.4 < 
PT < 0.6 GeV/c; (b) 1.4 < pr < 1.6 GeV/c 



ical pt bins from the data sample. The mass resolutions 
(a) for reconstructed A and A are typically about 3-4 
MeV/c^ based on a Gaussian fit to the peak. The back- 
ground beneath the A(A) peak is dominated by combina- 
toric pairs of charged particles. Decays of K'^ tt+tt" 
also contribute to the smooth background due to pions 
misidentified as protons. The yield is obtained from the 
invariant mass distribution in each transverse momen- 
tum (pt) and rapidity (y) bin, where the shape of the 
background near the A(A) peak is fit with a second order 
polynomial function. The signal to background ratios in 
the peak are greater than 6 for all the measured pt bins. 
Variations in the yield due to different fits for the back- 
ground have been included in the estimate of systematic 
errors. 

The raw yield for each pT-y bin was then corrected for 
finite detection efficiency. The overall correction factor 
was obtained from a Monte Carlo embedding procedure. 
A (A) particles were generated using a STAR GEANT 
package for detector simulation including all physical pro- 
cesses to take into account effects such as scattering and 
absorption in the detector support structures. A TPC 
Response Simulator (TRS) was used to generate TPC 
information at the pad level for all charged tracks inside 
the TPC from simulated A and A decays. The TRS in- 
cludes all TPC resolution effects ranging from electron 
transport in the gas to signal processing in the readout 
electronics. The simulated data were then embedded into 
real collision events from STAR detectors and were ana- 
lyzed with the standard STAR software package as real 
data. The correction factor for each px-y and multiplic- 
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ity bin was obtained from the number of A (A) particles 
reconstructed from the embedded sample, divided by the 
number generated. The combined acceptance and effi- 
ciency for A and A ranges from 0.8% to 5.8% as a function 
of pt for the most central collision sample. 
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FIG. 2: Transverse mass distributions of A (left) and A (right) 
at mid-rapidity (|j/|<0.5) for selected centrality bins. The 
dashed lines are Boltzmann fits. Note that multiplicative fac- 
tors have been applied to data from the two most central data 
sets for display. 
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FIG. 3: The mid-rapidity A (|2/|<0.5) transverse mo- 
mentum distribution from the top 5% most central colli- 
sions. For comparison the distributions for negative hadrons 
{(fN/{2TvpT)dpTd'n, |r?|<0.1) and anti-protons {\y\<0.1) for 
the similar centrality bin are included. Statistical errors are 
less than the size of the data points. 



Centrality 


0-5% 


5-10% 


10-20% 


20-35% 


35-75% 


dN/dy A 
(Boltz) A 


17.0±0.4 


13.0±0.3 


10.1±0.2 


5.9±0.2 


1.61±0.05 


12.0±0.3 


9.6±0.3 


7.4±0.2 


4.6±0.1 


1.26±0.04 


Tb a 
(MeV) A 


298±5 


304±6 


303±6 


289±6 


254±5 


312±6 


310±6 


305±6 


280±6 


258±5 


dN/dy A 
(exp fit) A 


17.4±0.4 


13.3±0.3 


10.4±0.2 


6.1±0.2 


1.66±0.05 


12.3±0.3 


9.8±0.3 


7.6±0.2 


4.7±0.1 


1.30±0.04 


Te a 
(MeV) A 


355±8 


364±9 


362±8 


343±8 


295±7 


374±9 


373±8 


366±8 


331±8 


301±7 



TABLE I: Fit parameters from Boltzmann and exponential 
fits of the TTiT spectra for A and A at mid-rapidity (|i/|<0.5). 
Only statistical errors are presented. The systematic errors 
on dN/dy and T are approximately 10%. 



The measured A spectra contain contributions from 
primordial A, S'^ decays, and feed-down from multiply- 
strange hyperons-notably S° and The primordial A 
and the decay products cannot be separated in our 
analysis and have been treated as primary A. We have 
investigated the feed-down from S hyperons by compar- 
ing systematically the distribution of the DCA to the 
primary interaction vertex for the reconstructed A from 
real data versus simulated primary A and secondary A 
from S decay. From this analysis, we estimate the con- 
tributions of feed-down from multiply-strange hyperons, 
mostly S and S decays, to be approximately (27 ± 6)% 
of the measured A and A yields respectively. 

Figure |2| presents the mr spectra (invariant distri- 
butions) of A and A for five selected centrality bins. 
Systematic errors from various methods of yield extrac- 
tion, reconstruction efficiencies and uncorrected sector- 
by-sector variations in the TPC performance are in- 
cluded. Both exponential (e-(™T-™)/'^^;) and Boltz- 
mann (mTe"'^'"^"'"^/^^ ) functions, as well as a Gaus- 
sian (e~PT'/(2'^ )) in have been used to fit the data, 
following the ansatz described in [ p5[ . 

All three functional forms were fit to the A and A spec- 
tra used for this analysis. The slope parameter obtained 
from the exponential fit is systematically higher than 
that for the Boltzmann by approximately 40-65 MeV. 
However, overall the integrated yields from all three fit 
functions are consistent within the statistical errors. The 
slope parameters and the dN/dy for the Boltzmann and 
exponential fits are presented in Table ^. The Boltzmann 
form was adopted in Figure |^ because it typically pro- 
vides a better ^-nd gives a reasonable description of 
the rriT spectra over the entire range of centrality and 
tranverse momenta investigated. 

Within the systematic error, the slope parameters mea- 
sured for the A and A tot distributions are the same. 
There is a systematic increase in the slope parameters 
from approximately 254 MeV for the least central (35- 
75%) to 312 MeV for the most central (0-5%) bin. Sim- 
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ilar behavior as a function of centrality is found for the 
p transverse mass distributions . Assuming the tem- 
perature at which particle interactions cease (the "freeze- 
out temperature pof) is constant independent of coUision 
centraHty, the increase in the slope parameter may be 
interpreted as an increase in the collective radial veloc- 
ity ^ Ip. 

A possible indication of hydrodynamic flow is the in- 
crease of the observed mean transverse momenta for vari- 
ous species with increasing particle mass. The transverse 
momentum distributions of negatively charged hadrons, 
as well as p and A are shown in Figure ||. The p and A 
Pt distributions are similar in shape in the region where 
they can be compared (below 1 GeV/c), even though the 
data sets cover somewhat different ranges in pT- Both 
distributions are qualitatively different from and much 
less steep than the corresponding h~ distribution, which 
is dominated by pions. Qualitatively similar behavior 
was observed in heavy ion collisions at the CERN Su- 
per Proton Synchrotron (SPS) In general, 
the slope parameters for all species are observed to be 
larger at RHIC, which is attributed to an increase in the 
collective radial velocity as a function of center-of-mass 
energy [ pO[. 

Figure^ indicates that at higher px (pt>1 GeV/c) 
the ratio of A to negative hadrons increases rapidly. The 
baryon to meson ratio at RHIC for pT>l GeV/c exceeds 
expectations from perturbative QCD inspired string frag- 
mentation models which were tuned to fit e'^e~ collision 
data and are the basis for modeling particle production 
in hadronic collisions as well |3^, For example, the 
A to h~ ratio is approximately 0.35 at pT of 2 GeV/c. 
Data from e+e^ collisions and calculations from string 
fragmentation models, however, indicate that although 
the baryon to meson ratio from quark and gluon frag- 
mentation increases as a function of Feynman x, the ra- 
tio never exceeds 0.2 As mentioned above, a natural 
explanation for the increase at high pT would be a large 
collective radial flow at RHIC Alternatively, 
it has also been suggested that the energy loss of high 
Pt partons could modify the baryon to meson ratio . 
To determine the exact dynamics that cause the rela- 
tive enhancement of baryons to mesons at high pT-, more 
experimental measurements over a larger p-p-range are 
needed. 

The average ratio of pt integrated yield A to A is 0.74± 
0.04 with no significant variation over the measured range 
of centrality. Given that there is a net excess of baryons 
at mid-rapidity, it is reasonable to conclude that there is 
more than one process contributing to A production and 
that significant baryon number from the colliding beams 
is transported to A hyperons at mid-rapidity. A question 
in that regard is why the shapes of m-r spectra for A and 
A are the same within errors. It has been suggested that 
significant rescattering of A and A during the evolution 
of the collision can lead to equilibration p9|, 54|. The 
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FIG. 4; A and A rapidity density as a function of negative 
hadron multiplicity at mid-rapidity. The open symbols in 
the figure are data points from an independent analysis of 
minimum bias data set using event-mixing techniques. Errors 
shown are statistical only. 



nature of baryon number transport itself requires further 
study. 

Figure ^ shows the dN/ dy of A and A from the Boltz- 
mann fit as a function of the h^ pseudo-rapidity den- 
sity 01 . At mid-rapidity the hyperon production is ap- 
proximately proportional to the primary h^ multiplic- 
ity in Au-|-Au collisions at RHIC. The dashed lines in 
the figure correspond to A = (0.054 ± 0.001)h^ and 
A = (0.040 ± 0.001)h- from a linear fit to the data. Sim- 
ilar centrality dependence of the lambda production was 
observed at the SPS energies. The A to h^ ratio at RHIC 
is much larger than that at the SPS while the A to h^ 
ratio is smaller at RHIC ||l^. This may be understood 
from the fact that at the SPS, most of the observed A 
hyperons carry baryon number transported from the col- 
liding nuclei through associated production, rescattering 
and fragmentation processes. In this case, the yield of 
A to h^ ratio is larger than that at RHIC due to the 
relatively high fraction of the A yield not resulting from 
pair production. Conversely, the increased importance of 
baryon pair production at RHIC energies must contribute 
to the observed increase of the A to h^ ratio relative to 
the SPS measurement. 

In conclusion, we have presented the first inclusive 
mid-rapidity (|y|<0.5) A and A spectra as a function of 
centrality from Au-|-Au collisions at the i/s^ = 130 
GeV energy. Salient features of the data include 1) large 
slope parameters of transverse mass spectra probably re- 
sulting from increased collective radial velocity at RHIC, 
2) similar shapes for A and A spectra despite the fact that 
a significant fraction of the A hyperons at mid-rapidity 
carry baryon number from the incoming nuclei while the 
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A hyperons are primarily pair produced, and 3) a signifi- 
cant increase in the A yield relative to negatively charged 
primary hadrons at moderate pt above 1 GeV/c which 
cannot be described by existing perturbative QCD in- 
spired string fragmentation models. The pt integrated 
rapidity densities of A and A are approximately propor- 
tional to the number of negative hadrons at mid-rapidity. 
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